IFUP-TH/2010-22 
SISSA 58/2010/EP 



Vortices in (2+l)d Conformal Fluids 



Jarah EVSLIN fl and Chethan KRISHNANt 



O ! ^ INFN Sezione di Pisa, Largo Pontecorvo 3, 

_ : Ed. C, 56127 Pisa, Italy 

^ ■ and 

Department of Physics, University of Pisa, 
Largo Pontecorvo 3, Ed. C, 56127 Pisa, Italy 

2 SISSA and INFN - Trieste Section, 
p i' Bonomea-265, 1-34136, Trieste, Italy 



Abstract 
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^ . We study isolated, stationary, axially symmetric vortex solutions in (2+l)-dimensional viscous 

O ' conformal fluids. The equations describing them can be brought to the form of three coupled 

flrst order ODEs for the radial and rotational velocities and the temperature. They have a 

rich space of solutions characterized by the radial energy and angular momentum fluxes. We 

K> ! do a detailed study of the phases in the one-parameter family of solutions with no energy flux. 

j^ . This parameter is the product of the asymptotic vorticity and temperature. When it is large, 

the radial fluid velocity reaches the speed of light at a flnite inner radius. When it is below a 

critical value, the velocity is everywhere bounded, but at the origin there is a discontinuity. We 

comment on turbulence, potential gravity duals, non-viscous limits and non-relativistic limits. 
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1 Introduction 

Incompressible fluids in two spatial dimensions allow vortex solutions with a rotational 
velocity profile that goes as 1/r. In this note we aim to generalize vortex solutions to the 
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relativistic case of conformal fluids with and without viscosity. We write down their equations 
in a tractable form, discuss their solution space, and study some of their phases in detail. 

What is the relativistic generalization of the 1/r rotational velocity? There are two natural 
guesses. If the velocity v continues to behave roughly as 1/r in the relativistic case, it will reach 
the speed of light at a finite inner radius. On the other hand if the relativistic velocity u = f 7 is 
roughly 1/r, then the velocity does not reach the speed of light at any finite radius. Interestingly, 
we find that both generalizations of the 1/r velocity profile are realized by relativistic vortices. 
When the product of the asymptotic vorticity and temperature exceeds a critical value the 
rotational velocity reaches the speed of light at a finite radius. Below this critical value the 
velocity remains bounded, although the inward radial velocity is nonzero and so discontinuous 
at the origin. 

1.1 Motivation 

Vortices are likely play an essential role in turbulence in any fluid. The picture of incompress- 
ible turbulence as composed of superpositions of vortices dates back to Richardson [1] , and is 
supported by evidence from analytical, numerical and experimental perspectives. The energy 
in such a turbulent fluid flows from modes corresponding to larger distance scales to smaller 
distance scales, in a phenomenon famously called "the cascade" . In 1941 Kolomogarov demon- 
strated the existence of certain scaling laws in this regime |1]. In incompressible (3+l)-d fluids, 
the conjectured regularity of the flows appears to be a consequence of vortex-stretching. For the 
incompressible Navier-Stokes equation in 3+1 dimensions, there exists a fully regular vortex so- 
lution that is stabilized by a radial velocity, called the Burgers vortex (and its generalizations). 
Vortex stretching is the phenomenon by which such a vortex with non-zero inward directed 
radial velocity can conserve volume because matter moves outward along the third dimension. 
These vortex tubes are believed to be the "sinews" of turbulence [3]. For compressible fluids 
in 3+1 dimensions, there is numerical evidence that vortices play a similar role (see e.g., [2]). 

The nature of (2+l)-dimensional vortices is quite different. The viscous dissipation may be 
arbitrarily small in 2+1 dimensions, and so one can work with the vortices of the non-dissipative 
(Eulerian) fluid, which have a divergent rotational velocity at the core. The statistical mechanics 
of these vortices (e.g.,|5]) is in broad agreement with experiments and numerical simulations^. 
In contrast, the (3+l)-d Burgers vortex depends on the viscosity crucially for its existence. 



"'^In "working with inviscid vortices, we ignore energy and enstrophy dissipation. The latter is a questionable 
assumption at high Reynolds numbers, but it has been argued that this assumption does not invalidate the 
results |6J. 

^The fact that there exists a regular vortex solution that is sustained by viscosity is undoubtedly of relevance 
to the existence and smoothness of the solutions to the Navier-Stokes equations. Burgers vortex has non- 



Another difference in (2+l)-d incompressible turbulence is the appearance of an inverse energy 
cascade [3 |H]: energy is transmitted to bigger and bigger length scales as vortices coalesce. 

Recently, the AdS/CFT correspondence [S] has opened up a new perspective on the old 
problem of fluids and turbulence. The claim is that the fluids that arise as the long-wavelength 
limits of (a large class of) conformal quantum field theories are expected to be described by 
gravity through the fluctuations of black objects in an asymptotically anti-de Sitter space, 
as described in Refs. pU| [TT| \T2\ [T3| |T3]. These references provide a map from reasonably 
generic fluid configurations to solutions of Einstein's equations up to some order in derivatives. 
Needless to say, if this implies that turbulence is as ubiquitous in gravity as in fluid mechanics, 
it would be difficult to overstate its relevance for example for cosmology. Indeed the chaotic 
evolution of gravity near a big bang singularity in the BKL ansatz is already well established 
[T5] . On the other hand it is unclear at the moment how useful this map is going to be for 
a direct understanding of turbulence, but the fact that Einstein's gravity can be dual to fluid 
dynamics is interesting in itself, and worthy of investigation. Therefore, because of the crucial 
importance of vortices in fluid dynamics, in this paper we initiate a study of vortex solutions 
in relativistic conformal fluids. 

One motivation for studying conformal fluids comes from the heavy ion collisions at RHIC 
and the LHC It seems possible that the state of matter that results from these collisions is 
best described by a strongly coupled phase of QCD plasma. This plasma is well-described by 
a turbulent [16] conformal fluid with a very small viscosity rj. AdS/CFT methods have been 
applied to study this system with very suggestive, but still only semi- quantitative, success [T7]. 
The onset of turbulence is expected to happen at high Reynolds numbers and for a relativistic 
fluid at temperature T, this is determined by 

sTL , , 

Re >1. 1.1 

Here s is the entropy density and r/ is the viscosity, and there are experimental and theoretical 
reasons to believe that r^/s is close to its conjectured lower bound, namely l/47r. For the 
RHIC plasma the relevant length scale L is the size of the gold nucleus, and the temperature 
is the QCD scale. Plugging in these numbers, one finds that Re ~ 150 ^ 1. This suggests at 
least naively, that one should expect that the plasma is turbulent, and a first step towards an 
understanding of this deeply non-linear regime is the study of vortices. In heavy ion collisions, 
the standard assumption is that one is dealing with a boost-invariant plasma in the longitudinal 
directions. This is the so-called "Bjorken flow" [18] and it makes the relevant fluid dynamics two 



vanishing velocity along the third axis. So contrary to what one might naively expect, the viscosity need not 
dissipate the flow to an equilibrium, as it would in compact domains with decaying boundary conditions. In 
(2+l)-dimensional fluids, the third dimension is of course unavailable, however vortices with velocities that blow 
up at the core are still a very useful approximation for understanding the statistical behavior of turbulence. 



(spatial) dimensional. This is one of the reasons why we focus on vortices in 2+1 dimensions. 
The (3+l)-dimensional case is also certainly of interest: in particular, it will be interesting 
to see if there are fully regular vortex solutions (like the Burgers vortex in incompressible 
fluids) in conformal fluids. See [19] for a different approach to turbulence via string theory and 
holography. 

Consequences of the potential existence of vortices in neutron stars have also received some 
attention [201 • However the vortices of interest contain (color) magnetic flux, and so the purely 
hydrodynamic vortices of the present note may be too primitive to shed light on their rich 
phenomenology. Vortices in neutron stars are generally treated in the superfluid approximation, 
however viscosity is necessary in a large region of the interior of a neutron star to avoid an 
r-mode instability that would be in conflict with the observed rotational velocities [21]. This 
provides one motivation for our interest in vortices in viscous fluids. 

Apart from their possible appearance in turbulent flows, vortices might also be relevant for 
an understanding of gravity solutions in AdS. Recently certain black funnel solutions which have 
horizons extending to the boundary of AdS were constructed [22l|23] using AdS C-metrics [21]. 
They were interpreted on the boundary as a fluid in equilibrium with a boundary black hole, 
possibly providing a holographic construction of Hartle-Hawking vacua for strongly coupled 
fleld theorieqj. A natural generalization would be to include rotation along an axis on the 
boundary. This corresponds to a spinning black funnel that has a horizon stretching all the 
way to the boundary. Of course, the spin has to be such that one does not run into trouble with 
superluminal rotation and similar pathologies. Constructing such solutions using generalized 
AdS C-metrics seems possible, see also [25]. The relation between fluids and gravity have 
generalizations away from AdS as well. The membrane paradigm for black holes suggests that 
black hole horizons can be treated in terms of fluid dynamics [261 [27] . There has been recent 
work on the construction of gravity solutions using the blackfold approach [2B] (generalizing 
previous work on the construction of exact vacuum solutions in higher dimensions [2S])- In AdS, 
the relation between the boundary fluid picture and the horizon fluid picture was connected 
via a Wilsonian RG flow in W, 



In fact a large class of the solutions that we flnd in this paper have a core of flnite radius 
where the fluid quantities diverge, while asymptotically they tend to those of a static black 
brane. This structure has similarities to what could be the dual of a spinning black funnel-like 
solution that might exist in AdS. We emphasize, however, that at this stage this is only an 
analogy. 

In this paper we survey vortices in a conformal relativistic fluid in flat (2+l)-dimensional 



■^The C-metric based approach of [35] does not work for constructing geometries where the boundary is flat, 
so these were done for geometries with hyperbohc spaces or compact spaces on the boundary. 



Minkowski spacetime. The appropriate generalization of this construction to other (2+1)- 
dimensional geometries which are hkely to have gravitational duals is left for the future. An 
explicit gravity dual will enable us to construct the boundary stress tensor directly and see how 
it compares with the fluid stress tensor. 

1.2 Summary of Results 

We numerically find axially symmetric vortex solutions which asymptote to a static fluid con- 
figuration. The latter condition means that we are studying isolated vortices, deferring the 
harder problem of vortex interactions. Correspondingly, from the perspective of the dual black 
objects in AdS, we expect the geometry to asymptote to a static AdS black brane. This should 
be dual to a static fluid configuration at constant temperature in the asymptotic region of the 
boundary. Indeed, in the black funnel cases that have been investigated, it is known that the 
fluid does tend to a static ideal fluid configuration in this asymptotic region |23]. In any event, 
the boundary conditions that we impose on our solutions are that the velocities should fall off 
at infinity and that the temperature should have a leading behavior that is constant. It should 
be noted that implicit in this assumption is the condition that the solutions are well-defined 
for large r. 

Our concern will mostly be viscous conformal fluids in this note, even though we also 
comment on the ideal limit. The equations become significantly more complicated when we 
include viscosity. One of the results of this paper is that these equations can be brought to the 
form of three coupled first order ordinary differential equations, one each for the temperature, 
the radial velocity and the rotational velocity (see Appendix B). In order to make statements 
about the asymptotic behavior, we will assume that the velocities and temperature can be 
expanded in inverse powers of r in the asymptotic region. This is a sufficient condition for 
asymptotic staticity, but is not strictly necessary. We find that there are many classes of such 
solutions that one can find, and that all of them that are asymptotically static have an inward 
radial velocity. In fact, this is not too surprising because in viscous fiuids rotation causes an 
outward radial energy-momentum transfer between the fiuid layers, which may be balanced by 
an inward radial mass fiow. If one allows no radial fiux of matter-energy or angular momentum 
in the vortex, then one finds that there are no asymptotically static solutions that are allowed 
by the fiuid equations. If on the other hand, one allows such fiuxes, then solutions are possible. 
We study one case in some detail, one where angular momentum fiows radially from infinity 
to a sink at the center of the vortex. It exhibits two phases, both of which have some form of 
singular behavior at the core, as in the incompressible case. This singularity can be at finite 
radius or at the center of the vortex, depending on the parameters. The rotational velocity 



divergecl at the core in the former phase. For vortices in the latter phase, the rotational velocity 
is finite everywhere, but the radial velocity is non-zero at the origin resulting in a discontinuity. 
The temperature diverges at the core in both phases. The radial component of the velocity 
(we always work in Landau gauge) is always nonzero and in fact is always inwards. This is 
analogous to the Burgers vortex, but of course without the extra dimension. 

Conformal fiuids at rest have only one scale, set by the temperature. The other quantities, 
such as the pressure, density and viscosity, are proportional to powers of the temperature. 
The constants of proportionality that we will use for numerics are those that arise from the 
correspondence to AdS black branes, however we will see that any other choice of coefficients 
simply leads to certain rescalings of our solutions. 

The hydrodynamics of a neutral relativistic fiuid is completely characterized by the con- 
servation of a stress tensor. In our axially-symmetric, stationary ansatz this consists of three 
ordinary differential equations, representing the local conservation of the components of the 
momentum 3- vector. In the case of the viscous, conformal fiuids of interest, these equations 
involve the temperature and its first derivative as well as the rotational and angular velocities 
and their first two derivatives. Therefore the solutions are characterized by five constants of 
integration. Fixing an overall scale and specializing to the case with zero radial energy fiux, 
one is left with 3 dimensionless parameters. We will find that only a 1-parameter family of 
these satisfies our asymptotically static condition. 

More explicitly, these conservation equations can be brought to the form of three first order 
equations: one for the radial velocity, one for the rotational velocity and one for the temperature. 
These equations explicitly involve two constants of integration, ci and C2, which we will relate 
to the radial energy and angular momentum fiux respectively. Clearly the counting of solutions 
of these equations must proceed as outlined above. In units where the velocity of light and the 
Newton's constant in the gravity dual are set to unity (cf. f l3.19tl3.20|) ). the vortex solutions 
are therefore determined by five dimensionful quantities, or four dimensionless oneqj. The case 
with no energy fiux from infinity corresponds to setting Ci = 0, and in this case, the asymptotic 
static constraint kills two of these. The remaining dimensionful scales can be taken to be the 
asymptotic temperature Tq and the vorticity u. When we speak of vorticity, we will mean the 
integrated nonrelativistic vorticity, which is the integral of the rotational velocity around the 
vortex on a contour far enough away that the velocity is much less than the speed of light. 
By rotational velocity we mean the velocity in the angular direction, which has dimensions of 
velocity and is equal to the angular velocity times the radius. 



''More precisely, the coordinate velocity diverges, while the proper velocity tends to the speed of light. 
^Note that the transport coefficients appearing in the equations are part of the description of the system, 
not the solutions. 



The only dimensionless quantity characterizing the (ci = 0) vortex is the product ooTq. The 
phase of the vortex and the velocity at its core may only depend on this combination. We will 
see that when a; To > 0.1, the rotational velocity reaches the speed of light at a finite radius, 
which we will call Tc. This radius therefore must be equal to u times a function of uTq. The 
(relativistic) radial velocity at Vc we will see is always inwards, and is equal to -\/3/2. The true 
radial velocity, which is the relativistic radial velocity divided by 7, therefore goes to zero at 
Tc- By velocity we will always mean the relativistic velocity unless specified otherwise. The 
rotational velocity diverges like (r — rc)~^/^. On the other hand when coTq < 0.1, the velocity 
never reaches the speed of light. The rotational velocity decreases to zero at the origin, while 
the inward radial velocity tends to a constant which varies from when coTq = to l/-\/2 when 
uTq ~ 0.1. The temperature at the core however diverges like 1/r. 

The divergences in the core are clearly unphysical. In the case of an incompressible fluid, 
in a physical flow they are always avoided. However the 1/r vortex solution is still useful 
because it serves as an idealization in terms of which the statistical properties of the flow 
can be described. Such divergences are avoided in real situations because it is impossible to 
accelerate a fluid particle to the velocity of light during a finite time t. Our solutions should 
instead be interpreted as the universal t — )■ 00 limits of vortices subjected to certain boundary 
conditions for a time t. At all finite times t the solutions are nonsingular. 

This argument also applies to conformal fluids. It should however be mentioned that in 
extreme regimes where the various derivatives are large compared to the temperature of the 
fluid, we should in principle include higher order transport coefficients in the fluid equations to 
determine the correct flow ^U\. But whether the flow goes to a regime where higher derivative 
effects become relevant depends on the initial conditions of the flow: a million dollars are 
available to those who solve a very closely related problem in a very closely related situation. 
This is not our concern here, we are interested in stationary configurations which can function 
as building blocks of a more complicated fiow. 

In the next section, as a simple warmup we write down the most general axially symmetric 
vortex solution for incompressible Navier-Stokes equations in 2+1 dimensions. The 1/r vortices 
that have been used in the statistical study of (2+l)-dimensional turbulence are a special class 
of these. Then we turn to a description of conformal fiuid dynamics and the equations that 
describe stationary vortices in them. In the remaining sections, we go on to study various 
kinds of vortices from various angles. We consider different phases of ideal and dissipative 
fiuids analytically determining their limiting behavior far from and near to the core and then 
numerically finding the entire solutions. While a quantitative understanding of the role played 
by vortices in (conformal) turbulence is lacking at present, we believe that an exploration of 
the space of these quasi-solitonic solutions is the place to start. 



2 Incompressible (2+l)-diniensional Fluids 

In the absence of forcing, incompressible flows in (2+l)-diniensions solve the Navier-Stokes 
equation 

^t^ + V ■ V^ = z/V^^, V-v = (2.1) 

where v = {vx-,Vy) is the velocity, ^ = dxVy — dyVx is the vorticity and v is the kineniatical 
viscosity. We will search for axially-symmetric vortex solutions of these equations in which the 
velocity asymptotes to a constant at radial infinity, corresponding to the ansatz 

K, Vy) = i-yfir) + xg{r), xf{r) + yg{r)). (2.2) 

In this section the functions /(r) and g{r) will stand for non-relativistic quantities, every- 
where else we will consider their appropriate relativistic generalizations. The incompressibility 
immediately forces a condition on the expansion function g: 

9ir) = ^, (2.3) 

where k is an integration constant. Substituting this back into the ansatz and using the vorticity 
equation of motion, we find a simple differential equation for / whose general solution is 

/(^) = ^ + ^2 + Csr'/^. (2.4) 

Equations (12. 4p and (12. 3p together comprise the general solution of the stationary unforced 
axially symmetric equations in 2+1 dimensionqj. Note that to obtain the actual radial and 
angular velocity, one must multiply g and / by r. 

It is evident that there may be solutions that fall-off at infinity if one sets C2 = 0, so long 
as the radial flow is either zero or inward (in the latter case k is negative and should satisfy 
I A; I > z/). In a non-relativistic theory, the function g directly captures the radial inward flux: we 
can integrate the radial velocity rg{r) around any circle and we get a flux ~ k. The statistical 
mechanics of the / ~ 1/r^ vortices is a useful model for two-dimensional turbulence [T, "5]. 



^There are some caveats to this statement. One is when there is no radial flow, fc = 0. Then the most general 
solution is of the form 

f{r)^^ + C2 + C3logr, (2.5) 

for any non-zero viscosity v . When there is no viscosity, but the radial velocity is non-zero (fc 7^ 0), then the 
general solution takes the form 

.fir) = ^ + C2. (2.6) 

When viscosity and radial velocity are both zero, any function /(r) is a solution. 



3 Equations of Motion for Vortices in Conformal Fluids 

In the rest of this note we will generalize the discussion of Sec. [2] to the case of stationary, 
relativistic, axially-symmetric vortices in a (2+l)-diniensional conformal fluid. Conformality 
implies that all of the transport coefficients of the fluid are equal to constants times powers 
of the temperature T. To determine these coefficients unequivocally, we need to have full 
control over the underlying CFT or the gravity dual. However we will see that changes in the 
values of these constants can easily be reabsorbed into a rescaling of the solutions. Thus, we 
lose no generality by choosing these constants to be those arising from the black brane ansatz 
considered in [13]. 



Because of axial symmetry, every object must transform covariantly under the SO (2) axial 
symmetry. The spatial velocity is an SO (2) vector. Therefore it must be of the form 

u"" = -yf{r) + xg{r), u^ = xf{r) + yg{r) (3.1) 

where u^ and u^ are the two spatial components of the relativistic 3- velocity, in other words, 
they are the true velocity multiplied by 



7(r) = Vl + (m")2 + (uyy = ^/l + r^iPir) + g^{r)). (3.2) 

We will use the word "velocity" for u. f{r) and g{r) are arbitrary functions of the radial 
coordinate, /(r) is the angular velocity, rf{r) is the rotational velocity and rg{r) is the radial 
velocity, m is a function of x and y, but the dependence will be left implicit. We will also 
sometimes omit the r dependence of /, g and 7. Indices will be raised using the metric r]^'' 
with signature (—1, 1, 1). In particular this yields the temporal component of the velocity 

u^ = 7(r) = -uo (3.3) 

and, being timelike 

M% = -1. (3.4) 

There is no canonical definition of the velocity -u of a fluid. We will use the velocity of 
the Landau frame, which is the frame in which the stress tensor consists of an isotropic non- 
dissipative term plus a dissipative term which is entirely spacelike. One may then obtain the 
spatial components in this frame by contracting with the projection 

p/.^ = r^^-^ + u^u" (3.5) 

/ r^iP + g^) {-yf + xg)j {xf + ygh 

i-yf + xg)-f 1+y^p + x^g^ -2xyfg xy{g^ - p) + (x^ - y^)fg 
\ ixf + yg)-f xy{g^- p) + {x^ -y'^)fg I + x^p + y^g^ + 2xyfg 



Being a projector, all eigenvalues of V are equal to zero or one. By construction the unique 
zero-eigenvector is the velocity u 

p/^^u^ = r]^'''u^ + u^'u^'u^ = n'^ - n^ = 0. (3.6) 

Also it is easy to see the V squares to itself. 

The dynamics of an uncharged relativistic fluid is simply described by the conservation of 
the stress tensor, which to first order in a derivative expansion is 

T^'^ = pu^'u" + PP^'^ - 77a^'^ - COV^"". (3.7) 

Naturally the derivative expansion is only reliable only as long as the derivatives are small 
compared to the temperature. 

P is the pressure. In a conformal theory the stress tensor must be traceless. We will define 
the pressure and density so that these first two terms, which are those of an ideal fluid, are 
traceless without the first derivative terms. Therefore 

p = 2P. (3.8) 

The third and fourth terms of (13. 7p are proportional to derivatives of the velocity, r] is the shear 
viscosity, and ( the bulk viscosity. 6 is the 3-divergence of the velocity, and a is a traceless 
symmetric tensor equal to the shear, which will be defined momentarily. The quantities p, P, 
rj and ( are all scalars, and so by axial symmetry can only depend on the radial coordinate 
r. The entire stress tensor must be traceless, and so we must also impose the tracelessness of 
the first order derivative terms. This implies that the bulk viscosity ( vanishes and the stress 
tensor may be simplified to 

T^"" = P{r]^"' + Su^u") - T](r^"'. (3.9) 

The matrix a is the shear, which is the spatial, traceless, symmetric part of the first deriva- 
tive of the velocity. To make it spatial, one need only multiply by the projector V, and to make 
it traceless, one subtracts the trace. One may subtract the trace before the projection, and it 
remains traceless because the trace of the terms projected out are of the form 

u^d.u^ = ^d,{u^u^) = ^9,(-l) = 0. (3.10) 

Therefore the matrix a is 

cr'"' = V^^'V'f'ido.Up + dpu^ - v^pdX) (3.11) 

where we have used the fact that there are two spatial dimensions in the normalization of the 
divergence term. 

10 



The component a** is an SO (2) scalar, the components a*'^ form a vector and a^'' is a 
symmetric 2-tensor. Therefore they may be written 

a" = a"(r), a*" = -A{r) + -B{r), a'^ = -A{r) - -B{r) 

7 7 7 7 

a"- = x^C{r) + y^D{r) + xyE{r), a"^ = {y^ - x')^ + xy{C{r) - D{r)) 

a^y = x^D{r) + y'^C{r)-xyE{r) (3.12) 

where as always 7 depends on r, but the dependence is left implicit, as will be the dependences 
of A, S, C, D and E. A few pencils later one can express these functions of r in terms of / 
and g 

a« = 2r^fgf + r\g^-f)g' (3.13) 

A = -r'f^f + g^)+r{l + r'g'){ff + gg') (3.14) 

B = -r\f + g^)fg-rgf{l + 2r^f + r^g^)+rfg\l + r^f) (3.15) 

C = -2f^ + i(l + rVy (3.16) 

r 

D = 2f^ + 2r/^/'-i(l + rW (3-17) 

r 

i? = 2f -2//--(l + rV)/'. (3.18) 

r 

In the stationary ansatz a** and i? do not affect the equations of motion. Notice that, as we 
have set the speed of light to unity, velocity is dimensionless and so / and g have the same 
dimensions as 1/r. The shear, a, is the derivative of a velocity and so it also has dimensions 
1/r, which implies that A and B have dimensions 1/r^ and C, D and E have dimensions 1/r^ . 

So far our discussion has been applicable to any conformal fluid. Any conformal fluid has 
a single scale, the temperature T. All other thermodynamic quantities may be expressed as 
monomials in T. In (2+l)-dimensions, tj is proportional to T^ and P to T^. 

T] = aT"^, P = bT^. (3.19) 

In fluid flow dual to a (perturbed) black brane in AdS, these constants of proportionality are 
fixed in terms of the Newton's constant G of the gravitational theory [13] . 

-^. 6-5^. (3.20) 

In fact, our numerical results are easily generalizable to the more general conformal fluid. In 
a generic conformal fluid, the deflnitions of t] and P in terms of T may be multiplied by some 

11 



numerical constants Ci and C2 respectively. We will see that if Ci = C2 then the equations 
of motion fl3.30|3.31ll3.32ll3.33p for the velocities are invariant, although both constants of inte- 
gration (ci and C2, cf. (I3.23p . (I3.28P ) must be multiplied by Ci. For general values of Ci and 
C2, one obtains the same equations but with the temperature T multiplied by C1/C2 and the 
constants of integration multiplied by Cf/C|. Therefore the velocity profiles that we find for 
general values of Ci, C2 and T can be straightforwardly translated to the general case. 

Putting together Eqs. dSH IXT^ IXM IXTBl IXT71 IXTHl KW. K^ we can express the stress 
tensor as a function of /, g and T. It will be convenient to decompose the stress tensor into 
SO {2) tensors, as we did for a 

T"" = x^C{r) + y^V{r) + xy£{r), T"^ = {y^ - x^)^ + xy{C{r) - V{r)) 

avy = x^Vir) + y'^Cir) - xy£{r) (3.21) 

and write the equations of motion directly in terms of the functions of the radius that appear 
in those tensors. The equations of motion are just the conservation of the stress tensor. As the 
configuration is stationary, the time derivatives vanish. 

Consider first the conservation of energy 

= (9^T^'* + dyTy' = 2A + rA. (3.22) 

This equation may be integrated to solve for A^ yielding 

A='j, (3.23) 

where c\ is a constant of integration. We will later interpret c\ as the radial energy flux, which 
falls into a sink at the core of the vortex when ci 7^ 0. This is easily seen by remembering the 
definition of the stress tensor as the flux of the /i-th component of momentum 3- vector through 
a surface at fixed v coordinate. The energy flux is computed as 

Flux = /" T^^'dY.^ ~ /" {T'^'dy - T'^dx) ~ I r^A d9 = 2-kci. (3.24) 

J T. JUr J r=const 

Because of the local conservation law we can move the integral to a circle of constant ra- 
dius, which we denote E^. We have changed from Cartesian to polar coordinates using x = 
rcosO, y = rsin^, while noticing that the surface element (which is one-dimensional in two 
spatial dimensions) goes as {dy,—dx). The final result arises from Stokes' theorem applied 
to (13. 22^ I3.23P if there is a source/sink at the origin with strength ~ Ci. We will mostly be 
interested in vortices with no net energy flux, so we will set ci = in most of this paper. 
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Next we will impose the conservation of momentum in the x direction 

= 4T- + 3„T" = x(3C + rC - V) + |(4£ + rS'). (3.26) 

This equation must be satisfied at every point in space, therefore the coefficients of x and y 
must vanish separately. This yields two equations 

P = 3C + rC (3.26) 

and 

A8 + r8' = Q. (3.27) 

The second equation may be integrated to solve for £, yielding 

^ = 5 (3.28) 

where C2 is a constant of integration which we will see is proportional to the product of the 
asymptotic viscosity and angular velocity. The equation for conservation of momentum in the y 
direction is related to (I3.25P by an SO (2) transformation, so it can give no new constraints on the 
S0(2)-invariant functions. An analogous computation to the one we did for the energy net-fiux 
can be done for momentum as well. In polar coordinates this will become angular momentum, 
with the resulting condition that C2 = should hold for the net angular momentum to be zero. 

The functions A^ B, C, V, and £ may be expressed in terms of the coefficients (13. 14113. 16f3.17f3.1^ 
by substituting their definition Eq. (l3.12p into Eq. (13. 7p . We will see that B does not appear in 
the equations of motion, the others are equal to 

T'' = (2 + 3r\f + g^))P-7]a'\ A = 3-fgP--A (3.29) 

7 

C = (^ + 3g^)P -vC V = {^ + 3f)P -r]D, £ = ~6fgP - r]E. 
Thus we may summarize by writing a complete set of equations of motion 

3^gP-!LA = % (3.30) 

C = (^ + 3g')P-r]C (3.31) 

{\ + 3f)P-r]D = 3C + rC' (3.32) 

6fgP + 7]E = -^ (3.33) 

where P and i] are given in Eq. (I3.19[ I3.20p in terms of the temperature and the Newton's 
constant in the gravitational theory. Eq. (I3.30p is the conservation of energy, Eqs. (I3.3ip 
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and fl3.32p are the conservation of momentum in the radial direction and Eq. fl3.33p is the 
conservation of angular momentum. The equations of motion are all proportional to the inverse 
Newton's constant. Therefore Newton's constant may be multiplied out of the equations of 
motion, and so does not affect the equations for the velocities and temperature although it 
provides a multiplicative constant in the viscosity and pressure. 

4 Ideal Conformal Vortices 

We will start by looking for vortex solutions in ideal fluids. The ideal fluid flow equations 
are obtained by setting rj = in the equations of the previous section. The equations that we 
need to solve take the form 

37sP-^ = 0, 6fgP + Sl = o (4.1) 

i{^ + ^'f)p]' + i{^-+3s')p] = {^,+3f-)P (4.2) 

An asymptotic expansion of the form (see discussion in section 15.11) 



/ = § + § + ••• (4-3) 



ly^Zi /ftj 



9 = ^ + ^ + ... (4.4) 



ly ^ (i^O 



yields the constraints 



P = P, + ^ + ^ + ^ + ... (4.5) 



Ci = 35-2^0, C2 = -6/2fi'2^0, (4.6) 



on /2 and (72 with no constraints on Pq. The higher order coefficients are determined recursively 
from these. This suggests that there are indeed asymptotically static vortex solutions in the 
non- viscous case. 

An interesting fact is that in the ideal case, these solutions can be analytically determined. 
The way to do this is to solve the first two equations (which are algebraic) explicitly and then 
solve the remaining first order ODE directly. The solutions are of the form 



where a is an integration constant and there is no correlation between the signs. Once / 
is known, g and P can be determined algebraically via the constraint equations. At first 
sight it might seem that these solutions do not allow asymptotically static solutions, but one 
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should remember that the integration constant can be imaginary as long as the final solution is 
real. Indeed it is possible to show that by choosing a appropriately, solutions that satisfy the 
asymptotic expansion can be found. 

When ci and C2 are both zero and we set g = 0,we can also see a different branch of solutions: 
any function /(r) is a solution with the temperature (and therefore pressure) determined by 

T{r) =Toexp ( / f{rfrdr\ (4.8) 

This function-worth of solutions is analogous to the solutions we found in the incompressible 
case, see footnote [61 

Now we turn to the discussion of vortex solutions in the presence of viscosity for generic 
values of ci and C2. 

5 Asymptotics and Phases of Viscous Conformal Vortices 

In this section we consider vortices in dissipative conformal fiuids by including the viscosity 
term and consider the full equations of motion that we found in Section |3l These equations 
consist of a first order equation (I3.33P , a second order equation (using (I3.3ip as definition of C 
in (I3.32P ) and a constraint that we will use to obtain the temperaturq^ (I3.30p . There are various 
kinds of non-linearities in these equations. An analytic solution is likely to be difficult to find, 
so we will try to find numerical solutions of this system. The trouble is that Mathematica fairs 
better with differential equations that are of the same order, rather than with a mixed system 
especially when there are constraints. So as a first step we will find a transformation that brings 
these equations to a set of three first order equations of the structurally very simple form 

/'(r) = F(/,^,T,r), (^'(r) = G(/,^,T,r), r(r) = i/(/,^,T,r), (5.1) 

for known (but somewhat complicated) functions F, G and H . The details and some explicit 
forms are given in Appendix B. 

We will seek solutions of these equations that are "asymptotically ideal" . That is, we want 
them to go to a static fiuid with constant non-zero temperature far away from the core. We 
discuss some features of the Ci 7^ 7^ C2 case first and then turn to a rather detailed analysis of 
the Ci = case in the following subsections. Finally we comment on the case where Ci = C2 = 
and discuss some subtleties involved in this case. Our aim is only to construct some solutions 
and not attempt an exhaustive scan of the solution-space. 



^The transport coefficients of a conformal fluid are powers of the temperature times a constant. 
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Figure 1: Plot of the temperature (curve non- vanishing at infinity), —g (the lowest curve) and 
/ for /2 = l,To = 0.5,(72 = —1- The horizontal axis captures the radial coordinate near the 
divergent core r^ ~ 6.998. 



When Ci = there is no energy fiux from infinity into a sink at the origin and when C2 = 
there is no such angular momentum fiux. Of course, in a realworld fiuid fiow singular sinks 
do not exist, but they are useful idealizations that capture dissipative features of the fiow. In 
the incompressible (non-relativistic) case, if there is no radial fiow there are no sinks: this is 
because material transfer is necessary for fiux of either kind to exist. The situation is more 
subtle in the relativistic case we consider. The condition for no fiux is not simply the vanishing 
of g: the precise conditions for no fiux are given by the vanishing of ci and C2 as we showed 
in the last section. In fact, we will see that asymptotically static vortices with both ci and C2 
zero are impossible in dissipative conformal fiuids (unlike in incompressible fiuids). 

5.1 A General Class of Vortices 

The most general isolated stationary vortex has both a fiux of energy and angular momentum 
from infinity into a singular sink in its core, corresponding to nonzero values of the constants 
of integration ci and C2. We plot a typical profile of the various vortex quantities as a function 
of the radius in Figure [H More exotic profiles are possible: one curious feature of these vortices 
is that there exist choices of parameters for which the vortex changes its direction of rotation 
at a finite radius, where it goes through zero rotational velocity. 

Asymptotic Expansion: We are interested only in asymptotically static (isolated) vortex 
solutions. This means that at radial infinity we want the velocity to tend to zero and the 
temperature to tend to a constant. It is useful (both for numerics and otherwise) to understand 
the mutually consistent behavior of the unknown functions at infinity that satisfy the equations 



16 



of motion. To do this we substitute the expansion 



/ = ^ + $ + - (5-2) 

9 = ;:2+^ + - (5-3) 

Ti To T, 

T = To + ^ + ^ + ^ + ... (5.4) 



into the ODE form of the equations presented in Appendix B and demand that they be zero at 
each order. There are no 0(^) terms for / and g because we want the velocities (rf and rg) 
to die down at infinity. Note that a solution of this form was strictly not necessar}o: but if we 
do find a solution, then that would be good enough for us (cf. Introduction). We find that the 
leading constraintqj on a consistent solution 

f -n _ n f - 9^2^ _ 9ciG 

h-0, g,-0, /2 - - (4^y^2 + i8,^c') ' ^'-A^' ^^-^^ 

A detailed analysis of the phases of the Ci 7^ 7^ C2 case would be interesting, but we will 
not undertake it here in this preliminary investigation. We demonstrate the existence of a rich 
phase structure by studying in detail the case where ci = 0, corresponding to no energy flow 
from infinity, in the remainder of this note. 

5.2 ci = 0: Large radius expansion 

It is straightforward to perform the asymptotic expansion as in the last subsection for the 
Ci = case as well, and fix the fall-offs. But instead we will proceed by direct inspection of 
the equations at large r. (Both methods lead to identical results of course.) First, inspired by 
the incompressible vortex presented earlier, we look for a rotational velocity rf that decreases 
as 1/r while g decreases faster at large r. This means that the angular velocity / decreases as 
l/r2 

/ = § (5-6) 

where /2 is a constant. The vorticity of a vortex is the integral of the curl of the velocity. In 
an ideal incompressible (2+l)-dimensional vortex, this is a delta function at the origin, and 
one may use Stokes' theorem to evaluate it as the integral of the rotational velocity along a 
fixed radius circle. We will adopt this terminology, and refer to the vorticity of a vortex as this 



^For example, exponential fall-offs at infinity are not captured by this form. 

^The subleading constraints can be thought of as equations determining higher order coefficients in terms of 
lower order ones. 
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integral, evaluated far enough from the core that the nonrelativistic approximation holds. The 
rotational velocity is then equal to /2/r and so the vorticity is 

u = 27r/2. (5.7) 

Now let us attempt to determine the leading order r-dependence of the energy conservation 
condition f l3.30p . P and 77 asymptote to constants due to asymptotic ideality. 7 tends to 1, as 
the velocity goes to zero. Thus the first term on the left hand side has the same r dependence 
as g, while the second term has the same r dependence as A, which is expressed in terms of / 
and g in Eq. fl3.14p . The dominant term at large r is 



u. 



A^rff'^-^ (5.8) 

which is of order 1/r^. On the other hand, the right side of Eq. (I3.30p is Ci/r^, which is of 
lower order in r than either A (which is of order 1/r^) or g (which falls off faster than 1/r^). 
In order for (13.301) to be satisfied, the order 1/r^ term on the right hand side must vanish, and 
so we learn that 

ci = (5.9) 

as desired. 

Now that Ci = 0, the two terms in Eq. (I3.30p must be equal 

37C/P = ^A. (5.10) 

As P, 7] and 7 all tend to constants at large r, g must have the same scaling as A 



92 

^4 ■ 



(5.11) 



The constant g2 is nonzero and in fact negative, as rj and 7 are nonzero and positive, and if the 
vorticity u is nonzero so is /2 by (15.70 and therefore A is strictly negative by Eq. (15. 8p . Thus 
we have found that the Landau frame velocity will always have an inward component, similarly 
to the ordinary velocity in a stretching Burgers vortex. 

From (15.1 op and (15.81) we see that this inward velocity is proportional to the viscosity and the 
rotational velocity squared and inversely proportional to the energy density. This is reasonable 
for the following reason. In the presence of shear viscosity, there is a net flow of angular 
momentum/energy from regions of smaller to greater tangential velocity: viscosity couples to 
velocity gradients. As the velocity decreases with radius, this means that energy-momentum 
flows from smaller to larger radii. In a stationary configuration, this outward energy flow must 
be canceled. It is canceled by an inward velocity. The inward energy flow is equal to the inward 



velocity times the energy density p = 2P. The content of (15.101) is precisely the balancing 
between these two terms. 



We may furthermore use (15.101) to find the ratio of the pressure and the viscosity 

P A 



(5.12) 



In Eq. (I3.19P we have already seen that this ratio may be expressed in terms of the temperature 

P AttT a 
V 3 372^ 

and so we may find the temperature in terms of the velocity 

A 



(5.14) 



This expression for T is exact, it is not the asymptotic value of T at large r. Such an exact 
result is possible because, once we have established that Ci vanishes at large r, it vanishes 
everywhere. We may use it together with (15. 8p to determine the asymptotic value Tq of the 
temperature in terms of the asymptotic velocities 

This is positive as g2 is negative. 

Note that the temperature is independent of Newton's constant G, it is written purely as 
a function of the velocities in the problem. However, as we used r] and P to arrive here, in 
conformal fluids without gravity duals this expression for T will generally be corrected by a 
multiplicative factor. 

Now that T has been expressed as a function of the velocities, it is easy to use Eqs. (I3.19[ 
I3.20p to express the other thermodynamic variables as functions of the velocities 

^" 1447rG7V' ^" 4327rG7V' ^^'^^^ 



One then finds the asymptotic values of 1] and P 



/4 _ / 6 
J2 p _ i2 

SQnGgi' ^~5AnGg'l 



where again Pq is positive because g2 is negative. 

Plugging the expressions (I5.16P back into the equations of motion (I3.30|3.31|3.32|3.33p one 
finds a system of differential equations that depends only the velocities, with no thermodynamic 
quantities. Moreover, by choosing t] and P as in Eq. (I5.16P we have satisfied the conservation 
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of energy condition f l3.30p . Eq. f l3.3ip is the definition of C as a function of /, g and their 
first derivatives, and thus may also be ehminated. Eq. f l3.32p depends on the first derivative 
of C and so is a second order ordinary differential equation in the velocities, while Eq. f l3.33p 
is a first order differential equation. Thus we are left with one ordinary first order differential 
equation and one ordinary second order differential equation for two functions /(r) and g{r). 
One then expects that solutions will be parametrized by C2 and 3 constants of integration. One 
of these is given by the vorticity u = 27r/2 and one by the asymptotic temperature Tq. The 
various phases of the solutions only depend on the dimensionless combination uTq. The fact 
that two of the constants of integration do not lead to parameters in our space of solutions is not 
surprising because typically we expect also growing modes to be present which are outside of 
our asymptotic fall-off ansatz. In our case in fact it seems also likely that the two missing modes 
might in fact be a pair of complex (conjugate) solutions for at least some regimes (of r and 
parameters) due to some suitably defined discriminant vanishing. We have found such complex 



solutions analytically for various special case^lj: they involve functions of the form a/C^ — r^ 
which are not well-defined for r > C for some integration constant C. We have not attempted 
a general analysis of such solutions because our interest is primarily in asymptotically ideal 
solutions. 

The constant C2 is easily evaluated in terms of the boundary conditions /2 and To and the 
conservation of angular momentum (13.331) . First one needs to find the asymptotic behavior of 
E from Eq. (I3.18p . It is dominated by the f'/r term 



r 
Substituting this into (I3.33P one finds 



f'. 
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C2 = -r\QfgP + TjE) ~ -r^TjE ~ -47^0/2 = -^^- (5-19) 

5.3 ci = 0: Vortices with bounded velocities 

As we report in Sec. [6l numerically we have found two different kinds of solutions in the 
cl = 7^ c2 case, characterized by different behaviors in the core of the vortex. One kind 
of solution, which occurs for ojTq < .1, is characterized by a bounded 3- velocity u. In this 
subsection we will analytically determine the leading order behavior of this solution at small 
values of r. 

The inward velocity rg tends to a constant value Wrad as r — ;■ while the rotational velocity 
rf tends to zero as a power of r 

f-fsr\ g-^, 6>-l. (5.20) 

r 



One of which we report in Subscc. 15.5 
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Recall that our differential equations for the velocity depend only on / and g, and so this is 
enough information to determine the velocities and also the temperature by Eq. f l5.14p . The 
pressure and viscosity may also be determined in terms of the dual Newton's constant G. 

One may now calculate the small r behavior of all of the functions that we have introduced. 
As the rotational velocity tends to zero, its contribution to 7 also tends to zero leaving 



7~V^ + "-d (5-21) 

and so not only does the relativistic inward velocity tend to a constant, but so does the true 
inward velocity M/7. Similarly the true rotational velocity tends to zero. 

The leading order limiting behaviors of the other functions follow similarly 



1 + ^ radJ'^rad --y ^ (^ + '"radJ'^rad ^ ^ "^rad j^ _ r,^_2 , x i X„,2 N f „5-2 



Ar.-'—^^^, C^- ' ^l^' "" , /^~^, i?~-2«<i + 5 + Kad)/.r 



^^21^' ^^-10^' ^-(-4M--2(^ + ^-L)ML)r^^^ (5.22) 

Notice that Eq. fl3.32p is satisfied, thus we need only check that angular momentum is conserved 
(13. 33 p . This equation implies that the leading order term in S is of order 1/r^ with a coefficient 
which is nonvanishing for a vortex with a; 7^ 0, in other words, for any vortex. Here we have 
found a leading order term of order r^~^. Therefore 6 may not be greater than 0, because there 
would be no order 1/r^ term. 

Thus we are left with the condition 

-1<5 <0. (5.23) 

Numerically we find that vortex solutions exist with all values of S in this range. Now the 
leading term in S must be order 1/r^, whereas the term in (I5.22p is of higher order. Therefore 
this higher order coefficient must be zero. Dividing through by —2fsu^^^ this implies 

2«L + '^^L + '5 = (5.24) 

and so 

«L = -^- (5.25) 

In particular the range of delta (I5.23P determines the range of Wrad and so 7 

- 1 < Mrad < 0, < 7 < V2 (5.26) 

where we have chosen a negative sign for qq in order to obtain a positive temperature. Thus 
the maximum true velocity M/7, which occurs for a vortex at the critical vorticity 

0.1 

Wc ~ 7^ (5.27) 
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is only 1/a/2 times the speed of light. On the other hand, the velocity at the core of the vortex 
may be relatively small if u is taken small enough and need not even be relativistic, although 
it is discontinuous at the origin. 

Although the velocity does not diverge, the thermodynamic quantities do diverge. Eq. f l5.14p 
and Eq. (15.161) imply that to leading order they diverge as 

2 3 

rp —Ura.d ^rad p ~^rad /r 2Q^ 

47rr ' '^1447rGr2' ^ 43271 Gr^' ^' ' 

Again these quantities are positive because Wrad is negative. These divergences are features of 
our stationary ansatz, and so are not expected to persist in vortices that have formed over a 
finite period of time. 

We have determined the leading order behavior of the solution at the origin in terms of 
the unknown variables Mrad and fs-, where /^ is increasingly irrelevant as it is the coefficient of 
the rotational velocity which tends to zero at the origin. Being dimensionless, u-^^ may only 
depend on the dimensionless combination of constants of integration wTq, we will write this 
fact as 

Mrad = Mrad(wTo). (5.29) 

We are unable to determine this function analytically, and so will find it numerically in Fig. [2] 
of Sec. O 

5.4 ci = 0: Vortices with a finite radius divergence 

For a given value of the asymptotic temperature Tq, if the vorticity u is greater than the critical 
vorticity Uc in (I5.27p . at a finite radius Tc the rotational velocity TcU is infinite, and so the true 
velocity is the speed of light. 

Numerically we have found that the behavior near r^. is universal. Rather than demon- 
strating that no other behavior is allowed, we will merely show that the behavior that we find 
satisfies the equations of motion. When r is slightly greater than Vc we find 

/ ~ r- I ' 3 ~ - — ■ (5.30) 

In particular, the angular velocity diverges. In this case not only the powers of the leading 
terms in the velocity but even the coefficients themselves are entirely determined by Vc- r^ has 
dimensions of length, like the vorticity u and so Tc/uj must be entirely determined by the only 
dimensionless combination of the constants of integration ojTq. Therefore 

Tc = uohiuTo) (5.31) 
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for some function /i, which is plotted in Fig. [3] of Sec. |6l 

Now that we have formulas for / and g, we know the limiting behaviors of the velocities and 
so, as in the previous subsection, may determine the limiting behaviors of the other functions. 
This time / dominates over g, and so 7 is dominated by the large rotational velocity 

7 ~ r J ~ -^^. 5.32 

In particular 7 diverges and so the true inward velocity rg/'-y tends to zero. The limiting true 
rotational velocity r//7 tends to the speed of light. 

The other functions are easily found to be 



(r-rc)2' rl{r-rcy r^ir-rcY' ^3/2(^ _ ^j3/2 

1 1 ^ 2-2 

C ^ , V ^ , £ :-r^ . 5.33 

3V37rGr2(r-r,)3' V^TiGrl{r - r,Y' TiGrl'\r - r,)^/^ 

Here the vanishing of the 0{{r — Vc)^'^^'^) term of 8 is expressed as above to show that, with 
the particular values of the constants in (15.301) . this leading term, which would be inconsistent 
with the equation of motion Eq. (13.331) . vanishes. In principle, to verify that (13.331) is satisfied 
with a nonvanishing value of C2, one should also calculate the order 1/r^ term, but we have 
simply verified numerically that this works. One can see from (I5.33P that the other equation 
of motion (I3.32p is satisfied. 

The thermodynamic quantities are also divergent at Tc- The temperature diverges as 

T ^^ (5.34) 

while the viscosity and pressure diverge as 

11 

"^ ^ 127rG(r-re)2' ^ 6v^7rG(r - rjs' ^^'^^^ 

The divergence of thermodynamical quantities in the infinitely boosted Landau rest frame of 
the liquid may not be particularly physically relevant. However the energy T** in the rest frame 
of the asymptotic fluid diverges even more strongly 

In fact there is nearly a cancellation leading to a lesser divergence, the viscosity term in the 
energy is minus one half of the pressure term. This semi-cancellation has some of the flavor of 
a Virial theorem. 
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Figure 2: When u < Uc, the velocity is everywhere bounded. The rotational velocity tends to 
zero at the origin, and the angular velocity diverges as r^. The leading order behavior may be 
entirely determined by Urad? which in turn is a function of the dimensionless combination uTq. 
In this figure we numerically determine this function. Notice that at the critical value of .1 the 
coefficient Wrad goes to —1, which we have seen is the maximum velocity allowed for a solution 
with this scaling. Beyond this critical value the rotational velocity diverges at a finite radius 
Tc, and the solutions are in the other phase. 



5.5 ci = = 02'. No asymptotically static vortices 

Finally we turn to the case Ci = C2 = 0, in which there still is no energy fiux but now also no 
angular momentum fiux. It is easy to see from the asymptotic behavior of the first order form 
of the equations in Appendix B that no asymptotically static solutions can exist (with non-zero 
temperature at infinity) in this case. But there is a slight subtlety to this statement. This is 
because in going from the original equations to the first order form, we have divided in places 
by g and so we are implicitly assuming that g y^ 0. This would have been a mere technicality 
had it not been for the fact that in the viscosity — )■ limit we did indeed find vortex solutions 
with g = in (14. 8p . This suggests that we look for an exact solution with g = 0. In fact, if one 
works with the original form of the equations of motion in the Ci = C2 = case and sets g = 
and looks for solutions, we can find the full analytic solution: 



T' 

— = fr, with / 



^/C^ 



(5.37) 



for some constant C. It is easy to see that this solution is not asymptotically static. 
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Figure 3: When u > Uc, the velocity reaches the speed of hght at a finite critical radius, rcrit- 
The dimensionless combination rcrit/to' is a universal function of the dimensionless combination 
uTq. This function is plotted above. As expected, the critical radius tends to zero as u 
approaches Uc, in accordance with the fact that the velocity is bounded when u < Uc- On the 
other hand at a fixed temperature in the high vorticity limit, the critical radius asymptotes to 
the vorticity. In particular for all solutions the critical radius is bounded by the vorticity. 
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Figure 4: This is the profile of the angular velocity / versus radius for the vortex characterized 
by the boundary conditions /2 = 10 and ujTq = 0.5. The rotational velocity reaches the speed 
of light at Tc ~ 1.208. At large radii the angular velocity scales as 1/r^, and near r^ it scales as 
(r — r^)"^/^. On the right we display the degree of /, defined as the exponent of (r — r^) in its 
leading term, or more precisely as (r — rc)f'/f. 
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Figure 5: This is the profile of the fractional contraction velocity —g, which is the true inward 
radial velocity times 7 divided by the radius, versus the radial coordinate. The vortex is 
characterized by the boundary conditions /2 = 10 and loTq = 0.5. While the rotational velocity 
reaches the speed of light at Vc ~ 1.208, the radial velocity tends to a constant. At large radii 
the fractional expansion velocity scales as 1/r^, corresponding to a radial velocity that scales as 
1/r^. On the right we display the degree of g, defined as the exponent of (r — r^) in its leading 
term, or more precisely as (r — rc)g' /g. 
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Figure 6: This is the profile of the degrees of the angular and radial velocity profile functions 
/ and g near the critical radius Vc ~ 1.208. These degrees are defined as (r — r^f'/f and 
(r — rc)g'/g. The vortex is characterized by the boundary conditions /2 = 10 and uTq = 0.5. 
One can see that the exponents of the (r — r^) scaling tend to the well-defined limits —0.5 and 
found in the asymptotic analysis. 
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Figure 7: This is the radial temperature profile of the vortex with uTq = 0.5 and /2 = 10. The 
temperature diverges at Tc, but at large r tends to a constant, with corrections of order 1/r^. 

6 Numerical Results 

Now we turn to a numerical investigation of the ci = vortices discussed above. Recall that 
in Subsecs. 15.31 and 15.41 we found two phases of solutions. The first phase is characterized by a 
finite inward velocity Wrad at the origin, while the rotational velocity tends to zero as r^+^. By 
dimensional analysis Wrad may be written as a function of the dimensionless combination uTq. 
In Fig. [2] we numerically find this function and we find that, as we have claimed, this phase is 
manifested when ojTq < 0.1. In particular, using Eq. fl5.25p . we see that as one approaches the 
critical vorticity Uc = O.I/Tq, 6 tends linearly to —1. This is to be expected, it implies that 
at the critical vorticity the rotational velocity no longer tends to zero. In fact, one reaches the 
limiting case Tc = of the second phase, in which the relativistic rotational velocity tends to 
infinity at r = r^. The critical radius of this phase, in turn, is plotted in Fig. [31 

We also display the profiles of some typical vortices in these two phases. First, we display 
graphical results for the numerical solution corresponding to /2 = 10 and uTq = 0.5. As uTq 
is five times the critical value, the velocity will reach the speed of light at a finite radius r^. 
In this case Vc ~ 1.208. The angular velocity / profile is displayed in Fig. HI We see that, in 
agreement with the asymptotic analysis in Sec. [5], at large r, / tends to f2/r'^, corresponding 



to a rotational velocity of /2/r. As r approaches r^ f diverges as l/-\/r — r^., corresponding to 
an angular velocity that approaches the speed of light. The radial velocity is seen in Fig. O 
Again, in agreement with the asymptotic analysis, the fractional radial velocity scales as 1/r^ 
at large r, corresponding to a radial velocity that scales as 1/r^. At large values of the radius 
this solution indeed becomes the usual nonrelativistic solution. Near the critical radius r^, g 
tends to a constant. This corresponds to a constant inward relativistic radial velocity. However, 
as 7 diverges, the true radial velocity tends to zero. In Fig. [6] we zoom in on the region near Tc, 
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Figure 8: This is the profile of angular velocity / versus radius for the vortex characterized 
by the boundary conditions /2 = 10 and uTq = 0.02, which is one fifth of the critical value. 
The rotational velocity vanishes at the origin as r"'^*^, and as always at large radii the angular 
velocity scales as 1/r^. 



seeing that the scaling exponents of / and g indeed tend to —0.5 and as claimed in Subsec. 15. 4[ 

The only dimensionful quantities that are needed to characterize the solutions are /2 and 
Tc, therefore if /2 is varied with ojTq fixed, all of these results must remain the same except that 
Tc will change such that /2/Tc is constant. Indeed this is the content of Eq. fl5.3ip . We have 
numerically verified that this is indeed the case. 

Finally we display a case in which the dimensionful parameter /2 is the same, but the 
dimensionless combination of asymptotic parameters is uTq = 0.02, which is one fifth of the 
critical value. Therefore the velocity remains everywhere bounded. In Fig.[H]we plot the angular 
velocity / versus radius, seeing again that at large r it scales as 1/r^, but that at at small r it 
scales as r~°'^^, indicating that the rotational velocity rf falls to zero as r*^'^* near the origin. 
The radial fractional velocity g is displayed in Fig. [91 Again asymptotically it shrinks as 1/r^. 
Now near the origin it diverges as 1/r, implying that the inward velocity rg tends to a constant, 
which is in this case 0.349. Thus 6 = —0.22 and the small r coefficient u^ad = 0.349, satisfying 
Eq. (15.251) . All of these results, being dimensionless, must not change if /2 is varied with the 
dimensionless combination uTq fixed. We have checked numerically that this is indeed the case. 
In Fig. [in] we see the radial temperature profile in vortex with uTq ~ .02. Consistently with 
the asymptotic analysis we see that the temperature tends to a constant at large r, with 1/r^ 
corrections while it blows up at small r. 
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Figure 9: This is the profile of the fractional contraction velocity —g versus the radius for the 
vortex characterized by the boundary conditions /2 = 10 and ujTq = 0.02. At the origin the 
inward velocity —rg tends to a constant value of 0.3485, while asymptotically it falls as 1/r^ 
as always. 
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Figure 10: This is the radial temperature profile of the vortex with uTq = 0.02 and /2 
The temperature diverges at the origin, but at large r tends to a constant. 
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7 Discussion 

Finding a vortex solution is clearly the most preliminary of steps in an approach to under- 
standing the dynamics of a fluid and in particular, turbulent flow. We have investigated 2+1 
dimensional vortices as a warm-up example: the intersting case of 3+1 dimensions should also 
be tractable. It would be interesting to see if one can find a solution that generalizes Burgers 
vortex to 3+1 conformal fluids. This might also be important in understanding real- world 
systems such as layers of neutron stars which are believed to be well-approximated by very 
low-viscosity conformal fluids. 

One of the questions that we have not investigated is that of the stability of the vortices. 
The fact that numerics is robust suggests that at the first order derivative expansion for the 
stress tensor that we are working with, the solutions are likely to be stable. At the core of 
the vortex where the quantities blow up, a regular solution might be possible by going beyond 
the first order terms. But as we have mentioned, stationary vortices serve as ingredients in a 
turbulent, time-dependent flow where singularities are avoided. In incompressible fluids, the 
full flow is expected to be understood purely using the second order (in the EOM, i.e., Navier- 
Stokes) terms for regular initial data. It will be interesting to determine for which regimes this 
approximation holds for conformal flows. 

An interesting aspect of the stability question is that due to the fluid-gravity correspondence, 
it could translate to questions of stability in black hole spacetimes. In fact, by taking a scaling 
limit of a conformal fluid flow, in [31], a map between the turbulent instabilities of Navier- Stokes 
at high Reynolds numbers and gravitational solutions was made. But there is a tautological 
nature to the fluid-gravity translation: for any flow satisfying certain general conditions one can 
write down a gravity solution order by order in derivatives. It is less clear at the moment how 
to put this to work and use one side of the duality to understand the long-standing questions in 
the other. Some related discussion on stability, higher derivative terms and regularity of flows 
can be found in the Introduction. Shock waves and boundary effects have also been discussed 
in this context in [32] l33] . 

7.1 Non-Relativistic Limit 

Asymptotically, the fluids we consider are dual to static black branes. The fact that the 
temperature is tending to a constant means that the (relativistic) density is also going to a 
constant. In fact it is easy to check that in the asymptotic region, our vortex solutions / and 
g also satisfy the incompressible Navier-Stokes equations up to corrections that die down at 
infinity. Using the incompressible fluid equations that we wrote down in section 2, for the 
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ansatz fl2.2p the incompressibility and Navier-Stokes equations take the form 

Incompressibihty = 2(7 + rg' = (7-1) 

Navier - Stokes = rg{2f + rf)' - -^ (r{2f + rf')] = 0. (7.2) 

r dr \ / 

It is easily checked that for the relativistic vortex solutions we found for the case Ci 7^ 7^ C2, 
in the asymptotic region 

Incompressibility ~ ^ ( " ) ' Javier — Stokes ~ Oi — ), (7.3) 

while for the case when ci = 7^ C2, both equations are order 1/r^. The velocities in both 
cases are of the order 1/r. Thus, our vortex solutions can be thought of as the UV (near-core) 
corrections to the incompressible vortex solutions. It will be interesting to see whether this is 
related to a certain scaling ( "non-relativistic" ) limit that was taken in [51] to reproduce the 
incompressible Navier-Stokes equation from conformal fluid dynamics. 

It is important to keep in mind that a constant density regime of the flow does not auto- 
matically mean that the relativistic nature of the underlying fluid can be ignored |37]. The 
reason for this is that in non-relativistic fluids, the pressure is always negligible compared to 
the energy density while this is not true for relativistic fluids. Our claim above is that the fluid 
velocities satisfy the incompressible equations of motion (up to some higher corrections in 1/r) 
as can be checked by direct computation. 

It should be noted that the question of what precisely constitutes a "non-relativistic" limit is 
somewhat ambiguous for relativistic fluids. In conventional fluid dynamics [35], the fluid particle 
is a black box with an equation of state and (local) thermodynamics, but with no microscopic 
structure. In particular, it is entirely consistent to view the thermodynamic variables like 
pressure, energy density and temperature as god-given quantities with no microscopic origin. 
This means that when one takes the slow-motion, i.e. non-relativistic, limit for fluid velocity 
(= velocity of a fluid particle), one has no canonical way of taking the non-relativistic limit of 
thermodynamical variables like pressure an temperature. To emphasize this, in an appendix we 
demonstrate a non-relativistic limit for barotropic fluids where the equation of state before and 
after the limit stays unchanged. This should be contrasted with the limit deflned in [34j (for 
the special case of conformal barotropic fluids), which changes the equation of state. In this 
latter limit, since the fluid becomes incompressible, the sound-propagation velocity becomes 
inflnite. 

Of course, it is possible to start with kinetic theory and some assumptions about the nature 
of the microscopic particles (eg., that they are monatomic and therefore there are no vibrational 
degrees of freedom) to derive the fluid dynamical equations from the Boltzmann transport 
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equatioio [36]. The various non-relativistic limits can then be systematically obtained via 
different scaling limits of the particle velocities as the speed of light is taken to infinity. But 
unfortunately, this strategy is of limited use for conformal fluids because kinetic theory assumes 
a particle description, which is necessarily weakly coupled, and not immediately useful for the 
strongly coupled (near-) conformal fluids at RHIC and elsewhere. But it might still be interesting 
and fairly straightforward to see such a construction for weakly coupled conformal fluids, in 
particular for theories like TV = 4 SYM which have marginal directions described by a tunable 
coupling. 
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Appendix A On Non-Relativist ic Limits of Barotropic Fluids 

In this appendix, we discuss taking non-relativistic limits of compressible relativistic fluids. 
This is of some relevance to our work because our vortices become non-relativistic in the 
asymptotic region. Also, we want to emphasize certain distinctions between what has sometimes 
been called a non-relativistic limit in the literature and a scaling limit. We will work with ideal 
relativistic fluid^ whose energy-momentum tensor we repeat here in a more convenient form: 

T''={^,+p)u'u'^ + ^,r- (A.l) 

The flow is described by its conservation. The dimensionality of the spacetime is not crucial, we 
will take it to be 3+1, with metric signature is { — , +, +, +}, so u^u^ = —1. We have included 
a factor of c^ in the deflnition of P because this is more natural in the non-relativistic limit. 
The precise deflnitions of P and p are meant so that they reproduce conventional deflnitions of 
pressure and temperature in the c — t- oo limit. We will deal only with barotropic fluids, which 
are fluids whose equation of state is of the form P = P{p). The conservation equation 

df.T^"' = (A. 2) 



^^ Typically this is done for non-relativistic set-ups, but a relativistic version is conceptually not very different. 
^^The non-relativistic limits of the kind we discuss here seem unexplored also in the case of dissipative fluids. 
It would be interesting to see how the ideas generalize. 
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gives us four equations. The fluid is described by P,p and the four components u^. Together 
with the identity u'^u^ = —1, and the equation of state, this is enough to fully determine the 
fluid§ 

Ideal non-relativistic fluids are described by the continuity equation and the Euler equation. 
In the barotropic fluid case, the equations take the form 



dp 

— + V ■ (p v) = (Continuity equation) (A. 3) 

(9v VP 

— — h (v • V)v = . (Euler's equation) (A. 4) 

These give rise to four equations and we also have the equation of state P = P{p)- Counting the 
three components of v, the pressure P and density p we again have equal number of unknowns 
and equations. Even though we use the same symbols for the pressure and density in relativistic 
and non-relativistic fluids, they are a-priori logically distinct entities. The precise relation will 
be clarifled when we discuss the precise way in which the non-relativistic limit is taken. 

Our aim is to formalize the transition from flA.2l) to flA.3IIA.4]) . One approach for doing this 
for the specific case of conformal fluids was taken in ^34j to relate conformal relativistic fluids 
to the incompressible Navier-Stokes equation. We will follow a somewhat different approach. 
It seems quite likely that what we present here was well-known to the ancients, but we will 
present it here because our perspective might be slightly different. 

For the stress tensor given in flA.ip . the fluid-flow equations flA.2l) can be transformed to 
the compact form 

Vp+(p + ^]dX = 0, (A.5) 



V^d^P ( P 






Vv!" = 0. (A.6) 



We have deflned 



V^" = g"" + u^u" , V = u^'d^. (A.7) 



There is some redundancy in this description, as can be seen by contracting the second equation 
with Ufj, and noticing that it yields an identity. 

To make the non-relativistic limit intuitive, we will flrst write the relativistic equations in 
non-relativistic notation. To this end, we can identify u'^ = 7(1, v/c), where 7 = 1/a/1 — v'^/c^, 



^^In general a fluid equation of state can include other thermodynamical state variables: it can include for 
instance, T. This addition of an extra variable will mean that we will need an extra equation to determine the 
fluid flow completely. This extra equation is usually a thermodynamical relation between the state variables. 
See [33]. 
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which differs from the 7 defined in the text. This leads to the following useful explicit expres- 
sions: 



ax = ^Dv' + ^V ■ V, V = ^D, (A.8) 



---2^^^^ + ^--' ^-? 



where D is the convective derivative dt + (v. V) that appears in the Euler equation. Using these 
expressions and suitable linear combinations of (]A.5|) and (]A.6p . after some work, we can write 
down the relativistic fluid equations as 

d,p + V- (pv) + i^,[P+-,)Dv' + ^^ = 0, (A.IO) 

•^v , „, 1 / VP V a^P \ ,, ^^, 

The advantage of this particular form is that the structure of these two equations is an imme- 
diate generalization of the non-relativistic continuity flA.Sp and Euler f lA.4p equations. Now it 



is easy to see what we need to do in order to reproduce those equations from the relativistic 
ones. The appropriate limit is one where v, P and p and their derivatives are held fixed, while 
we send c — )■ 00. The equations immediately reduce to ( 1A.3P and (lA.4p . 

We see that the non-relativistic limit allows us to keep P and p independent, so the speed 
of sound Cs = ydP/dp is held fixed as c — > 00. This is essential, since otherwise the fluid 
would always end up being incompressible in the non-relativistic limit [3^. Note also that the 
P that appears in the incompressible non-relativistic fluid equations (see eg., equation (17) in 
|34] ) is in fact P/p in our (non-relativistic) notation. The non-relativistic limit taken in [31] 
for conformal relativistic fluids, P = p/3, was a different scaling limit, which involved retaining 
other combinations (that include c) of quantities fixed, while sending c — )■ 00. We think that 
the non-relativistic limit presented here is another very natural scaling limit. 

An interesting possibility is to see how the symmetries of the system change when we take 
our non-relativistic limit. Note that as a special case of the general thing we did above, we 
can consider the conformal case {p = p/3). The result is actually a non-relativistic fluid with 
p = p/3. This is not incompressible. It would be interesting to see what the conformal invariance 
of the parent fluid becomes in this limit. For discussions on this issue, see 
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Appendix B The First Order System 

In this appendix we will re-write equations (13.30113 .3^ as a simple first order system. In 
principle, we can write P and 77 in f l3.30p in terms of the T using f l3.19tl3l20|i^i . and then solve the 
constraint for T explicitly in terms of the fluid velocities and use it in the remaining differential 
equations. We have done this, and found that the results match in this approach with the 
approach presented below in some specific examples, as a check of our numerics. But in general 
this direct approach is cumbersome because the constraint is cubic and the resulting solution 
is very complicated, and the system is of mixed differential order. 

Instead we will work with the linear (in derivatives) system where the expressions are in 
fact small enough to be presentable. First, we write P and rj in terms of T as before. Now, 
observe that using the definition of A, (I3.30p can be re-arranged to the form 

g' = X,if,gJ',T,r). (B.l) 

(We won't write explicit forms until the very last step.). Note that (I3.3ip is just the definition 
of C. Combining (I3.3ip and (I3.32p . we find that the second order differential equation that one 
ends up with is of the form 

g" = X,{f,gJ',g',T,T',r) (B.2) 

This expression is linear in T'. The third and final equation is (13.330 . which can be written in 
the form 

f' = F{f,g,T,r). (B.3) 

Taking a derivative of (IB.ip we find that the derivative of Xi should equal the right hand side 
of (lB.2p . This gives us a solution for /" of the form 

f = X,{f,gJ',g',TX,r). (B.4) 

The right hand side of this in turn can be equated to the derivative of the right hand side of 
(lB.3p . The resulting equation can be solved for T', and we find 

T' = X,if,gJ',g',T,r). (B.5) 

At this stage, we can take the independent equations to be 

f' = F{f,g,T,r), g' = X,if,g, f,T,r), T' = X,{f,g, f',g',T,r). (B.6) 



^'^For a generic conformal fluid, the other thermodynamic variables are determined by temperature only up 
to proportionality constants. We will leave these constants arbitrary. 
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Using the first of these to substitute for /' in the RHS of the second, and then using the 
resulting equation as well in the third, we can bring it to the final form 

f' = F{f,g,T,r), g' = G{f,g,T,r), T' = H{f,g,T,r), (B.7) 

which is the promised first order system. The explicit forms of F, G and H are 
C2 + 2hr^fT'^ - 2r^fgT'^{hg - 3aT) 



G 



H 



2r3 {p + 5^2^2) r^2 

(ca/ + 2ci7 - ^br^Pg^T^ - Gar^gT^ - Qar^g^T^) 
2br^gT^ + 2br^g^T^ ' 

4c?7 + AciC2r^f + 02^7^^ - Uacir^gT'-^ - 2Aab-fr^ g^T^ + 

-Uar^g^T^ (ci + 2b-fT^) + 2c2f (2ci + 2cir^g^ - 3a-fr^gT'^) + 
+P (7 (C2^ + 4ci V2) - l2acir^gT^ - 2Aab-fr^ g^T^) 
2Aabr^g^ (1 + r'^g'^) T^ 



Note that we have kept the transport coefficients general in these expressions. To go to the 
black brane values, which is what we use in the numerics, we can use f l3.20p . In the above 
equations 7 = ^/l + r^^P + r^ as before. 
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